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We consider the hydrodynamics of relativistic conformal field theories at finite temperature and 
its slow motions limit, where it reduces to the incompressible Navier- Stokes equations. The symme- 
tries of the equations and their solutions are analyzed. We construct exact solutions with finite time 
singularities of one-dimensional relativistic conformal hydrodynamic motions, and use them to gen- 
erate multi-dimensional solutions via special conformal transformations. These solutions, however, 
are shown to have no non-trivial slow motions limit. A simple non-equilibrium steady state in the 
form of a shock solution is constructed, and its inner structure is analyzed. We demonstrate that 
the derivation of the gravitational dual description of conformal hydrodynamics is analogous to the 
QQ ] derivation of hydrodynamics equations from the Boltzmann equation. The shock solution is shown 

, to correspond to a domain-wall solution in gravity. We show that the solutions to the non-relativistic 

' incompressible Navier-Stokes equations play a special role in the construction of global solutions to 

! gravity. 
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\Q . I. INTRODUCTION 

The hydrodynamics of relativistic conformal field theories has attracted much attention recently, largely in view 
^ of the AdS/CFT correspondence between gravitational theories on asymptotically Anti-de-Sitter (AdS) spaces and 
O ^ CFTs [1] (for a review see 0)- Hydrodynamics gives a universal description of the large-time dynamics of the theory: 
O ] starting from an arbitrary initial state of the CFT, within the correlation time Tcor the system approaches the state of 
i-^^' local thermal equilibrium. At i 3> Tcor the evolution is mainly the evolution of the parameters of the local equilibrium 
described by hydrodynamics. On the basis of the AdS/CFT correspondence one then expects that the large-time 
dynamics of gravity can be obtained as a dual description of the CFT hydrodynamics [1, 0] ■ 
^ ' A complete (compressible) hydrodynamics is described by five fields: the three velocity components, the temperature 
^sD , and the particle density [EBB- ^ CFT there is no locally conserved charge corresponding to the particle density 
' @, Q. As a result, conformal hydrodynamics is described by only four fields: the three velocity components and the 
temperature. 

Hydrodynamics applies under the condition that the correlation length of the fluid Icor is much smaller than 
the characteristic scale L of variations of the macroscopic fields. In order to characterize this, one introduces the 
. dimensionless Knudsen number 

OO 
O 



Kn = Icor/L . (1) 

Since the only dimensionfuU parameter is the characteristic temperature of the fluid T, one has by dimensional 
analysis, 

Icor = {hc/kBT)G(X) , (2) 

where A denotes all the dimensionless parameters of the CFT. The function G(A) characterizes the CFT. 
The stress-energy tensor of the CFT obeys 

d^Tt"" = 0, T/^ = . (3) 

The equations of relativistic hydrodynamics are determined by the constitutive relation expressing T'^" in terms of 
the temperature T(x) and the four- velocity field u^^{x) satisfying m^m^ = —1. Here and T determine the local 
thermal equilibrium distribution of the fluid. The constitutive relation has the form of a series in the small parameter 
Kn < 1, 

CSO 

T^-'{x)=Y,Tr{x), Tr-{Kn)\ (4) 
/=o 

where Tj^'^ {x) is determined by the local values of and T and their derivatives of a finite order. The smallness of 
Ti'^ arises because it involves either the l—th derivative of or T or because it contains the corresponding power 
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of a lower-order derivative. Keeping only the first term in the series gives ideal hydrodynamics, within which the 
entropy is conserved and the entropy density per unit volume a obeys a conservation law d^{<Tu^) — 0. Dissipative 
hydrodynamics arises when one keeps the first two terms in the series. 

The ideal hydrodynamics approximation for T'"' does not contain the spatial derivatives of the fields. The stress- 
energy tensor reads (up to a multiplicative constant) 



(5) 



where r]^i, = diag[—, +, +, +]. 

The dissipative hydrodynamics is obtained by keeping the I = 1 term in the series in Eq. ^ . In the Landau frame 
d, Q the stress-energy tensor reads (up to a multiplicative constant) 



(6) 



where ct^j/ obeys a^^u^ — and is given by 
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(7) 



The dissipative hydrodynamics of a CFT is determined by only one kinetic coefficient - the shear viscosity rj. The bulk 
viscosity ^ vanishes for the CFT, while the absence of the particle number conservation and the use of the Landau 
frame allow to avoid the use of heat conductivity, which is not an independent coefficient here [8] . Dimensional analysis 
dictates that r/ = F{X)T^, where F{X) is a function characterizing the CFT and we again omit a multiplicative constant 
in Ti"'. 

The hydrodynamics of relativistic conformal field theories is intrinsically relativistic as is the microscopic dynamics. 
However, it has been shown in [10], that the limit of non-relativistic macroscopic motions of a CFT hydrodynamics 
leads to the non-relativistic incom pres sible Euler and Navier-Stokes equations for ideal and dissipative hydrodynamics 
of the CFT, respectively (see also [llj). For ideal hydrodynamics the implication follows by noting that the equations 
([3]) written in terms of v defined by u'' = (7, jv/c) and a variable P equal to InT up to an additive constant, take 
the form 



2c2 



dP 
'dt 



^ + (. • V).. 



-Xv-V)P 



1-^ 
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(8) 
(9) 



Now if we consider the solutions for which v remains finite in the limit c — s- 00, then the limiting field v obeys 



dv 
'dt 



+ {v ■ V)v = -VP - 



V • t> 1 I dP 



(10) 



The first equation above implies that VP must be finite at c ^ 00. However it is still possible that P contains 
a function of time leading to a finite spatially constant divergence of velocity. The general form of non-relativistic 
dynamics is therefore 



dv 

— + (v ■S/)v ^ -\/P - a(t)v, \/-v = -3a{t). 
ot 



(11) 



If we now impose the condition that v remains finite at large distances (which is valid in most physical situations), 
then we must require that a{t) = 0. We conclude that non-relativistic, finite at infinity, motions of the CFT obey the 
incompressible Euler equations 



dv 

— + (?;• V)v ^ -VP, V • t; = 



(12) 



Analogous considerations hold for the viscous hydrodynamics where the limiting field v obeys the incompressible 
Navier-Stokes equations 



dv 
'di 



+ {v ■V)v ^ -VP + i^V^v, V-t> = 



(13) 
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The kinematic viscosity v is given by 

V = nc'^F{X)/4:kBTo . (14) 
To is the main component of the temperature, which is approximately constant in the considered limit, 



V c 



(15) 



see [Toj . Thus, a relativistic conformal field theory contains the incompressible Euler and Navier-Stokes equations 
inside it. 

The solutions to the non-relativistic incompressible Navier-Stokes equations have special importance in light of the 
AdS/CFT correspondence: they allow to construct an explicit approximate solutions to the five-dimensional Einstein 
equations with negative cosmological constant 

+ 4g„,„ = 0, R=-20, (16) 



where we use RAdS = 1- If 'W and P solve Eqs. p3)) and obey the condition Kn <C 1, then the metric go defined by 

{h^kB = 1) 

igo)mndy"'dy'' = ~2Uf,{x")dx>'dr -f TT^T^{x")r-^u^{x")u^{x")dx>'dx'' + r'^ri^^dx^'dx" , (17) 
( I vie \ P 

has a value of Rmn + '^gmn, which is small both in Kn and v/c (see the details in section IVl) . 

Consider a typical solution of Eq. (llSp . which is turbulent (the majority of flows in nature are turbulent). There are 
two important scales in the solution: the outer scale of turbulence Lq, dictated by the boundary or initial conditions, 
and the Kolmogorov scale 

Ik - {v^/ef'^ « Lo , (19) 

where e is the energy dissipation per unit volume At the scale Ik the flow becomes smooth (here we make a 
rough estimate neglecting intermittency), so the Knudsen number obeys Kn ~ I cor /Ik- Using e ~ v^/Lq where Vc is 
characteristic value of the fluctuating component of the flow velocity, one may reexpress Kn in terms of Lq and Vc as 

Kn^ilcor/Loy^HvJcf/' . (20) 

Here we used the fact that generally F{X) ^ G(A), see [l3, [13 • sees that one can realize Kn ^ 1 by dialing 

the outer scale of turbulence Lq, and the characteristic fluctuating component of the velocity Vc- The possibility of 
having a globally defined hydrodynamic flow with a uniformly small Knudsen number is specific for non-relativistic 
flows, and it normally does not exist for relativistic (compressible) flows, see details in sections ITVl and FVl 

In this paper we consider the hydrodynamics of relativistic conformal field theories at finite temperature and its slow 
motions limit, where it reduces to the incompressible Navier-Stokes equations. The symmetries of the equations and 
their solutions are analyzed. We construct exact solutions with finite time singularities of one-dimensional relativistic 
conformal hydrodynamic motions, and use them to generate multi-dimensional solutions via special conformal trans- 
formations. These solutions, however, are shown to have no non-trivial slow motions limit. A simple non-equilibrium 
steady state in the form of a shock solution is constructed, and its inner structure is analyzed. We demonstrate 
that the derivation of the gravitational dual description of conformal hydrodynamics is analogous to the derivation of 
hydrodynamics equations from the Boltzmann equation. The shock solution is shown to correspond to a domain-wall 
solution in gravity. We show that the solutions to the non-relativistic incompressible Navier-Stokes equations play a 
special role in the construction of global solutions to gravity. 

The paper is organized as follows. In the section 2 we consider the symmetries of relativistic conformal hydrody- 
namics and their relation to the symmetries of the incompressible Euler and Navier-Stokes equations in the limit of 
slow motions. In section 3 we construct exact solutions of ideal CFT hydrodynamics. We analyse one-dimensional 
macroscopic motions, introduce Rieniann variables and provide a simple wave solution that becomes singular in a fi- 
nite time. We use the solution in order to generate multi-dimensional solutions via special conformal transformations. 
We show, however, that these solutions have no non-trivial slow motions limit. In section 4 we discuss weak solutions 
of ideal hydrodynamics defined at all times. We construct a stationary shock solution as an example of a global 
weak solution. We describe its inner structure due to viscosity, its conformal transformation and its non-relativistic 
limit. In section 5 we argue that the derivation of the gravitational dual description of conformal hydrodynamics 
is analogous to the derivation of hydrodynamics equations from the Boltzmann equation. We show that the shock 
solution corresponds to a domain-wall solution in gravity. 
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II. SYMMETRIES OF CFT HYDRODYNAMICS 



The conformal group is the symmetry group of the CFT equations ([3]), and it maps sohitions into solutions. The 
group includes the Poincare group, a dilatation and four special conformal transformations. The Poincare group 
reduces in the limit of slow motions to the Galilean group, including Galilean transformations and space and time 
translations, which is a symmetry group of the Navier-Stokes equations. 

Consider next the dilatation that acts as: 

Ax'^, ^ m''(Ax"), r(x") ^ X-^T{\x") . (21) 

The dilatation is a symmetry of the CFT equations ([3]), and it remains a symmetry for ideal hydrodynamics in the 
limit of non-relativistic macroscopic motions, i.e. a symmetry of the Euler equations. It acts as 

t^Xt, x'^\x\ v'^v'. (22) 

For viscous hydrodynamics in the limit of non-relativistic macroscopic motions, one has the dimensionfuU kinematic 
viscosity that breaks this symmetry, and indeed, ([22]) is not a symmetry of the Navier-Stokes equations. 
The Euler equations have a bigger symmetry Q , which is generated by (|22p together with 

t nt, x' x\ ^ ^-^v' . (23) 

However, the symmetry (|23|) does not follow from the symmetry group of relativistic conformal hydrodynamics in the 
limit of slow motions, cf. below. A combination of the two symmetries (|22|) and ([23]) with the same transformation 
parameter X = fi = a yields 

t ah, x' ax\ a-^v^ , (24) 

which is a symmetry of the non-relativistic viscous equations. In view of the above discussion, we conclude that also 
this symmetry does not follow from the symmetry group of relativistic conformal hydrodynamics in the limit of slow 
motions. 



A. Special Conformal symmetries of relativistic motions 

The special conformal transformations of the space-time coordinates x^ are given by 



^ ' l + 2b-x + b^x^ ^ ' 

where is a constant four-vector. They satisfy the identity 

.-(a,a>")(a.a>.) . . (26) 

A main distinguishing property of the CFT hydrodynamics is that the special conformal transformations allow to 
generate new solutions out of a given one. Consider a stress-energy tensor T^^{x) of the conformal field theory, i.e. 
that satisfies equations ([3]). Then, special conformal transformations produce a four-parametric family 



where each Tfj^^{x,b) satisfies equations ([3]). 

The vanishing of the trace T^{x, b) follows immediately from Eq. ([^^ . It is also straightforward to check that the 
equation dyT^'^{x, 6) = is satisfied, where the derivative is taken with respect to x. Thus, conformal transformations 
generate from any solution of ([3]) a four-parameter family of solutions. Consider now the relativistic conformal 
hydrodynamics, where the stress-energy tensor is defined by the four fields: and T. The analysis above implies 
that special conformal transformations of any solution u^, T of the conformal hydrodynamic equations generate a 
four-parameter family of solutions T{x,b), u^(a;,5), where 

Uf,{x, b) = {l + 2b-x + 52x2)(a^$")u„ mx, b)] , (28) 
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One can verify that Uf^{x, h) has the correct normahzation Uf^(x, b)u^{x, b)r]^'^ = —1 and the above transformations of u 
and T lead to transformation (|27p both for the ideal hydrodynamics given by Eq. ([5]) and for viscous hydrodynamics 
T^" given by Eq. ([6]). Note that the velocity transformation ([28|) can be rewritten as 

uf'ix, b) = [$(a:, b)] + 2u" [<S>{x, b)] 

baXf^ - b^'Xa + 2baX^'b ■ X - b'^XaX^ - x'^bab^^ 

l + 2b-x + 62^2 • ^^^> 



B. Special conformal symmetries of non-relativistic motions 



In order to analyze the non-relativistic limit of special conformal transformations of the hydrodynamics it is con- 
venient to consider the special conformal transformation of v defined by Ua = {—j,"fv/c). Using Vi = —cui/uo, one 
finds that Eq. (pS)) implies the transformation of v 



i{x,b) = 



{\ + 2b-x + b^x^)vi mx, b)] + 2vj [$(a;, b)] [Vx^ - b^x' + 2b>Xib ■ x - b^x^x, - x'^Vh 



l + 2b-x 



2 [2b^xab ■ X + 
2c [b°x^ 



bh 



x^bl] 



2vj b)] [bixo - baxi + 2b^xab ■ x - b'^x^XQ - x^h>bo] /c 



hx° + 2b°Xib ■ X - b'^x°x^ 



'b%,: 



l + 2b-x + b'^x^ + 2 [2b°xob ■ x + b'^xl + x^fog] _ 2vj [<^{x, b)] [W xq - b^xi + 2&Jxo6 • x - b'^x^x^ - x'^Wb^] jc 



.(31) 



Consider solutions of relativistic conformal hydrodynamics, where v remains finite in the limit c — > oo, and thus obeys 
the non-relativistic incompressible Euler (Navier-Stokes) equations. Then, if we employ in the above transformation 
such parameters b that v(x^ b) is also finite at c — > oo, we find a symmetry of the Euler equation. First it is necessary 
that the argument ^(x,b) is finite in the limit c — > oo. This implies 6* = a^/c^ and b^ = a/c where a* and a have 
either finite or zero limit at c — s- oo. Thus, 



t 



t 



14-2 



I -at' (l-ai)2 ■ 

The expression pip gives the transformation of v in the non-relativistic limit 



v{x, t, a,a) = V 



t 



at'' 



1-at' (l-at)2 



at — ax 



1 



at 



If we take a = (or equivalently assume that a scales as some negative power of c) , while is finite at c 
we get the following symmetry of the Euler equation 

v{x, t,a) = V [x — at^, t] 



(32) 



(33) 



oo, then 



2at 



(34) 



This symmetry, discussed also in |11| . corresponds to adding a uniform gravitational field to the Euler equation. 
Indeed, if v{x,t) satisfies the incompressible Euler equation then v{x,t,a) satisfies dtv + {v ■ V)v = — Vp + 2a. 
However 2a = V(2a • r) can be reabsorbed in the pressure, so that v{x,t,a) is also a solution of the incompressible 
Euler equation. Thus, this symmetry corresponds to a combination of a passage to a uniformly accelerating frame and 
the possibility of reabsorbing a (uniform) gravitational field into the pressure (known in fluid mechanics see e.g. pp. 
16 — 17 in [iSl). The above consideration also works for the Navier-Stokes equation. 

If we choose a^ — (equivalent to assuming that decreases as a negative power of c), while a remaining finite at 
c ^ oo we get 



v{x, t,a) = V 



t 



X 



1-at' (l-ai)2 



2ax 
1 - at' 



(35) 



This transformation describes an addition of a total expansion (from the center at a; = 0) or a total contraction (to 
the center at a; = 0) to the given flow. However, this transformation does not preserve the condition of finite velocity 
at large distances and as a result is not a symmetry of the incompressible Euler equation. In particular, out of an 
incompressible velocity it produces a field with V • v{x, t, a) = 6a/ (1 — at). 

Although the transformation (|35p breaks the finiteness of velocity at infinity, and does not describe a symmetry of the 
incompressible Euler equation, it is a symmetry of the more general equation pip . Note, also that this transformation 
is different from the special conformal transformation that completes the Schrodinger symmetry group. The latter 
acts as 



t ■ 



t 

1-at' 



1-at 



(36) 
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III. EXACT SOLUTIONS OF IDEAL CFT HYDRODYNAMICS 
A. One-dimensional motions and finite time singularities in ideal CFT hydrodynamics 

Equations ^ can be written as Q 

I?e = -^<9.«^ Pu^ = -a''e+^^^, (37) 

where ^ = InT and T) = u'^da. The first equation fohows from the second by multiplication with and the use of 
-u^-u^ = -1. 

Let us consider one-dimensional hydrodynamic motions, that is motions that depend only on one spatial coordinate, 
say X, and have = = 0. It should be stressed that the microscopic dynamics of the system is still (3 + 
1)— dimensional, and it is only the macroscopic motion that is one-dimensional. For example, one could think of 
pushing the CFT fluid through a pipe where x is the coordinate along the pipe. It is convenient to parameterize the 
velocity by u'^ — cosh(/) and = sinh0. Substituting u'^ = (cosh(/), sinh(/), 0, 0), we find that the /j, = and yU = 1 
components of the second of Eqs. (|57|) give 

d( smh2(j) d(f> cosh 2(j) — 2 dcj) 

^ ' c TT ' (38) 



dt 3 dt 3 dx 

cosh 2(l) + 2d(f) sinh 2<p 



dx 3c dt i dx' 

respectively. This system has a more symmetric form when written as an evolution equation for ^ and 

9^ _ csinh(20) dS, c dcj) 

'dt ~ ~ [cosh(2</.) + 2]'dx ~ [cosh(20) + 2]'dx' 
d(f> csinh(20) dcj) 3c d£, 

'dt " ^ [cosh(2(?!)) + 2]'d^ ~ [cosh(20) -I- 2] ft^ ' 



(39) 

(40) 
(41) 



It is interesting to notice the particular form that the above equations take, when considered up to quadratic order 
in (j). Introducing w'/c — 20/3 and p = we find 

These equations describe an isothermal one-dimensional flow of gas with temperature c^/3. They also describe one- 
dimensional isentropic, i.e constant entropy density, motion of a barotropic gas with velocity v' and pressure p = pc^ /3. 
In general, a barotropic gas with isentropic motion obeys p = Ap'^ where A is a constant and 7 > 1 is the polytropic 
exponent of the gas, which is equal to the ratio of the heat capacities at constant pressure and volume, respectively. 
Here 7=1, which corresponds to gas in the limit of a large number of internal degrees of freedom of the molecules. 
It is notable that a motion of such a gas admits a gravity dual interpretation, as we will discuss in the last section. 
The system of equations is a 2 x 2 system of quasi- linear PDE. It can be recast using Riemann variables, 



which are conserved along characteristic directions in space . Here the result of the procedure can be readily seen 
directly from the equations. Multiplying the second of Eqs. (|40p - ([^T1) with and adding and subtracting it from 

the first equation we get 



d£, , c [sinh2(/)± ^/3] d£, _^ 1 



90 c [sinh 20 ± ^/3] d(j> 
dt cosh 20 + 2 dx 



0. (44) 



dt cosh 20 2 dx V3 
These equations relate derivatives of and ^ along the two families x±{t) of characteristic directions, defined by 

dx± c [sinh 20 ± 
dt ~ cosh 20-1-2 



(45) 



Thus, two families of characteristics emanate from each point of the fluid. They describe non-linear propagation 
of sound to the right (x+)and to the left (a;_). In the limit of linear sound (small 0), the equations are simply 
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x± = ±c/\/3 where c/a/3 is the speed of sound. Correspondingly, the Riemann variables have a remarkably simple 
form 



V3' 



(46) 



and the equations governing one-dimensional hydrodynamic motions of the CFT fluid can be written in the form 



d sinh[V3(r+ - r_)] + ^/3 d 

dt cosh[\/3(r+ - r_)] + 2 dx 

d sinh[V3(r+ - r_)] - ^3 d 

dt cosh[%/3(r+ - r_)] + 2 dx 



= 



The above system admits special solutions, called simple waves, where one of Riemann variables is a constant, 
example, for r_ — const^ we find a closed equation on r+ of the form 

dtr+ + f{r+)d.^r+ = . 

Introducing w = f{r-\-) we find that w obeys the Hopf equation 



dw 
'dt 



dw 
ox 



(47) 

(48) 
For 

(49) 
(50) 



The equation describes formation of discontinuities (shocks) in w in a finite time. The shocks form when the spatial 
derivative of the initial condition on w is negative somewhere. Indeed, introducing a = dxW we find that a obeys 



dt 



a[xit),t] = -a^[x{t),t]^ 



dx 



w[x{t),t] 



(51) 



This equation predicts finite time blowup of a, that occurs at x{t) emanating from the initial point where a has the 
(negative) minimum. It is easy to see that the blowup corresponds to the formation of a discontinuity in w. One can 
use the implicit form of the solution of Eq. ((50|) , 

w{x, t) = wti[x — w{x, t)t] , wq = w{x, t — 0). (52) 

Qualitatively, Eq. (|50p describes velocity field of free particles moving on a line where the discontinuity forms when 
faster particles catch up with the slower particles in front of them. 

The formal reason for the formation of discontinuities in a simple wave is the intersection of characteristics. The 
field w in Eq. ([5D|) is conserved along the characteristic lines x{t) defined in Eq. ([5T|) . In general, nothing prevents 
the lines from intersecting in a finite time T. At t > T the solution obtained by imposing the field conservation along 
the characteristics would become multi-valued, which is signaled by the divergence of the spatial derivative of w at 
the intersection point. The same property of the finite-time formation of infinite derivatives of the fields holds for the 
solutions of the system (|T7|) -(|i5 1) and thus, also of the original system Indeed, the equations imply that 

r± are conserved along the corresponding characteristics in [x, t) plane. Generically, these characteristics intersect in 
a finite time, and at the intersection time a field derivative diverges. Thus, we conclude that ideal hydrodynamics 
equations in one dimension lead to a finite time singularity of the derivatives. 



B. Special conformal symmetry generation of higher-dimensional solutions 

As shown in the previous section, one may solve rather generally the one-dimensional dynamics of the fluid in a CFT 
hydrodynamics. We can generate from the one-dimensional solutions, higher dimensional ones by applying special 
conformal transformation as in Eqs. (P5)) - I|29p . With the help of Eq. ([501) . oii*^ finds that a solution of one-dimensional 



8 



hydrodynamics leads to a family of solutions: 

u°{x, b) = cosh (I) [^{x, b)] 



1 + 2 



{bx° - xb" 



1 + 26 • x + 



2sinh(?!)[$(a;,6)] 
1 + 26 • a; + b'^x'^ 



bix" - 6°xi + 26ia;°6 • x - b^xix^ 



-xV 



, u-^(a;, 6) = sinh</i [$(a:, 6)] 



26V6-x-62(xi)2-a;2(6i)2 



+2 



1 + 26 • X + 62^2 

box^ — b^xo + 2box^b ■ x — b^xox^ — x^bgb^ 



2cosli(/)[$(a;,6)] 



u^{x, 6) 



1 + 26 • x + 62^2 
2 



1 + 26 • a; + 622;2 



cosh0 6)] 60X2 - 62X0 



+2boX2b ■ X — b^XQX2 — 2^^6062 J + sinh (j) [^{x, 6)] 

X (61X2 — 62X1 + 261X26 • X — 6^x1X2 — x^6i62) 
e(x, 6) e [^{x, b)] + In [1 + 26 • X + 62x2] , 



(53) 



(54) 



where b^ = (6",^) and x^ = (x°,a;), the expression for is obtained from the expression for by interchange of 
indices, and and ^ are obtained from the one-dimensional hydrodynamics. 



C. The non-relativistic limit 

The relativistic one-dimensional solution does not have a non-trivial slow motions limit since the incompressibility 
condition implies that the velocity is independent of the coordinate x. This is in accord with the presence of a 
factor of c in front of the space derivative in Eqs. ([171) - (|15)) . We have seen that we can generate from the relativistic 
onc-dimcnsional solution, higher dimensional ones by applying special conformal transformations. It is still natural to 
ask whether these higher-dimensional solutions have a non-trivial slow motions limit, i.e. whether they also solutions 
of the incompressible Navier-Stokes equations. We will show in the following that this is not the case, and also the 
higher-dimensional solutions do not have a non-trivial slow motions limit. Of course, this is not really a surprise, since 
as we noted above, the one-dimensional solution with which we started in order to generate the higher-dimensional 
ones, has a trivial limit. 

For simplicity we set the time coordinate to zero. One can write the general structure of the velocity using 
equations ((54)) as 



w,(r) = (1 + 26 • a; + 62a;2)0($)V,$ , 



where 



b^x 



1 + 26 • a; + 62a;2 



The incompressibility condition implies then that 



Vi [(1 -t- 26 • a; + 62a;2)v,$) 
(l-H26-a; + 62a;2)(V$)2 



(55) 



(56) 



(57) 



Since the LHS of ((57|) is a function of <&, we need to require that the RHS is also a function of <i> in order to be able 
to solve of (j>. However, a straightforward calculation shows that the RHS of (|57p is not a function of $, for $ finite 
in the slow motions limit. Thus, the incompressibility equation cannot be solved, and there is no non-trivial slow 
motion limit to the higher-dimensional solution of the relativistic conformal hydrodynamics. 



IV. SHOCK SOLUTIONS 



A. Shocks - global solutions in ideal hydrodynamics 



We have seen that the relativistic ideal conformal hydrodynamics evolution in one dimension produces discontinuities 
in a finite time. The mechanism for the generation of the singularity is the intersection of characteristics, and it holds 
in higher-dimensional motions as well. A simple concrete example is the multi-dimensional generalization of the Hopf 
equation ([50)1 provided by dtW + wVw = 0. The matrix aij = djWi obeys the equation da/dt = —a^ , that predicts 
a finite time blow up of cr. The finite-time generation of discontinuities is a general property of ideal hydrodynamics. 

A natural question that arises is how should we define the hydrodynamics evolution beyond the singularity time? 
One way is to include the viscous contribution to the stress tensor, and consider Eqs. ^ and ([6]) instead of Eqs. ([3]) 
and ([5]). According to the ideal hydrodynamics equations, the viscous contribution grows in time and becomes infinite 
at the time of formation of the discontinuity [cf. Eq. ([7])]. Thus, the viscous contribution is important near the 
singularity. It is conjectured (cf. [1,01), that this contribution regularizes the singularity, makes all the fields finite 
and provides a well-defined evolution within the frame of differential equations. 

A different approach for dealing with the singularity is to continue the solutions of the ideal hydrodynamics beyond 
the singularity time, by passing from classical (differentiable) to weak solutions of the equations, where weak solutions 
are solutions with discontinuities. This formulation is very useful when considering high Reynolds number flows, for 
which the viscosity is effectively small Q. In such flows the viscosity is important only near the points, where a 
discontinuity would form in the absence of viscosity. The effect of viscosity is to smooth the discontinuity, while in 
the other regimes of the flow the viscosity effect is small. 

An appropriate framework is to apply an ideal hydrodynamics evolution everywhere, except at the hypersurfaces 
where a discontinuity forms. At the discontinuity one requires that the normal components of fluxes of locally 
conserved energy and momentum densities remain continuous. This requirement is necessary in order for the equations 
to describe the local conservation laws, as the hydrodynamics should do. The corresponding constraints at the 
discontinuity are called the Rankine-Hugoniot conditions. 

The main advantage of the obtained weak solutions - the shocks - is that they provide within the framework of 
ideal hydrodynamics, global and valid solutions at all times. Let us construct a simple example of such a solution. 
The simplest solution of the system of Eqs. ^ and ([5]) is a boosted equilibrium state, where T and u'^ are constants. 
We can saw two equilibrium solutions along a plane, and form a simple but non-trivial solution - a shock wave. T and 
are piecewise constant and they have a single discontinuity along a plane x = x{t). In the frame that moves with 
the shock where x{t) = (the considered shock solutions move at the constant speed), the fluxes of the conserved 
quantities must be continuous through the shock. This leads to the Rankine-Hugoniot conditions 



(58) 



where the brackets stand for the value of the jump in the function through the shock. In contrast, the entropy flux 
does not have to be continuous - rather the entropy grows in the shock. 
Using = (7, 71;) (here we set c = 1), we get from Eq. 



0, 



rp4 



1 -i;2 



, 



(59) 



and 



^T^VyVx 



1 



0, 



l-z;2 



0. 



(60) 



If Vx is non-zero the equations imply that the velocity components parallel to the shock plane are continuous across 
the shock. For convenience, we can set these components to zero, and we are left with 



1 



= 0, 



1 



= . 



(61) 

AaT'^ (where 



Note, that if we use the general conditions for the shock [6^ with p = aT^ and heat function i 
a = const), then we get the same conditions. 

If we denote by Ti , T2 and vi , V2 the values of the temperature and the velocity component at two sides of the 
shock, then the conditions read 
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These equations define a two-parameter family of shock solutions of hydrodynamics. The two parameters correspond 
to the values of the fluxes of the energy and the x— component of the momentum through the shock. Evidently, 
Ti — T2 and vi = V2 solve the above equations. In order to analyze the conditions for a non-trivial solution, we 
assume that the fields T2, W2 at a; < are given (one can assume t;2 > as the solution with V2 < can be obtained 
by the sign reversal). Equating the expressions for T2/T1, we find a quadratic equation for vi that has two solutions, 
vi = and the trivial solution vi = V2. Thus, a non-trivial solution with vi =^ < 1 exists only for V2 above the 

threshold value of 1/3 (recall that the speed of sound is l/\/3). In the frame for which the fluid is at rest at —00, the 
shock always moves at a speed larger than 1/3. Note that the difference vi — V2 can be arbitrarily small. 

Fixing the temperature from Eq. (p^ . we find that non-trivial shock solutions at V2 > 1/3 are described by the 
relations 

-1 = ^, T,4 = T,4|f^. (63) 
3u2 3(1 - 1^2) 

It is remarkable that the Rankine-Hugoniot conditions for the GET have such explicit solutions, cf. the usual situation 
0. Note that vi > V2 for 1/3 < W2 < 1/V3 while vi < V2 for V2 > l/\/3 where l/-\/3 is the speed of sound. Eor V2 
close to the speed of sound we have vi « V2 while at the exact equality V2 = l/VS the solution trivializes. Einally vi, 
like V2, belongs to the range 1/3 < wi < 1 as should be due to the symmetry in the indices interchange. 

The above shock solutions are stationary non-equilibrium states of the fiuid with finite dissipation in the shock. 
The stationarity is sustained by the fact that the solution is non- vanishing at infinity, and the presence of finite fluxes 
in the system. 



B. Special conformal transformation of shocks 



Shocks are quite simple solutions of ideal hydrodynamics: they are time- independent, one-dimensional and piecewise 
constant. As a result they allow a clear insight into what kind of a change the special conformal transformations 
produce in the solution. Using Eqs. ((3T|) and ([29]) wc find that shocks produce the following weak solution of the 
hydrodynamics 

(l + 2b - X + b'^x'^)v!k]Sii + 2v!k] \b^Xi — biX^ + 2b^Xib ■ x — b'^x^Xi — x'^b^bA 

Vi[x,b)- 



l + 2b-x + b^x^ + 2 [2b"xob ■ x + b^x^ + ^^b^] - 2v[k] [b^xo - b^x^ + 2b^xob ■ x - b^x^xo - x^b^bo] /c 
2c [b°x, - b,x° + 2b°x.,b ■ X ~ b'^x°x, - x^b%i\ 
'l + 2b-x + b^x^ + 2 [2bOxob ■ x + b'^xl + x'^bl] - 2v[k] [b^x^ - box^ + 2b^xob ■ x - b'^x^xo - x^b^bo] jc 



(64) 



where [fc] = 1, 2 and vy^^ T^fc] are some constants obeying Eqs. (|62p . The equation of discontinuity hyperplane (where 
one should switch /c = 1 to fc = 2 in the above formula) is given by x^ -I- b^x • a: = 0. One observes that the special 
conformal transformations allow to generate quite complicated motions out of even relatively simple motions. 



C. Shocks in viscous CFT hydrodynamics 

The ideal hydrodynamics shock solution constructed above is clearly not uniformly valid in space. It contains a 
discontinuity region, where derivatives are formally infinite indicating that one must account for the viscous, and 
maybe also higher order derivatives contributions to T^^ . The ideal hydrodynamics shock describes correctly the 
exact solution far from the region of fast variations of the fields, and it arises when the formal limit of zero viscosity is 
taken in the exact solution. Here we provide the exact solution of the viscous hydrodynamics, and demonstrate how 
the one-dimensional shock considered in the first subsection arises in the limit. 

The stationary one-dimensional solutions of the hydrodynamics equations satisfy 

d^T^"" = aoT^° + a^^T^^ = 0, aoT^° = . (66) 

These equations imply that 



T"^ = Ci, T^^ = C2 



(67) 
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where Ci and C2 are the constant fluxes of energy and momentum, respectively. 

As aheady mentioned above, the states in question are non-equiHbrium steady states sustained by the presence of 
finite fluxes in the system. In an ideal hydrodynamics, the non-trivial solutions of the above equations with different 
values of velocity and temperature at x = — c» and x — +00, are obtained by introducing a discontinuity at x = 0. 
In contrast, in the viscous hydrodynamics the different asymptotic values at a; = —00 and x = +00 are connected 
smoothly and the discontinuity at a; = changes to a smooth transition from the values of the fields at a; = —00 to 
those at a; = +00. Note, that at large |x| the derivatives are small and ideal hydrodynamics is a good approximation. 
The regions that are infiuenced by the viscosity are at the core of the shock. 

In order to take into account viscosity, we consider cr^^ in Eqs. ©-([T]). One notes that the time-derivatives 
appearing in ct^j^ can be substituted by their ideal hydrodynamics expressions. Indeed, within hydrodynamics the 
time derivatives of fields are given by a series in the spatial gradients of the fields. The second derivative terms in ct^^ 
are higher order terms in the expansion compared to the considered order involving the first derivatives fl5:|] . The use 
of Eqs. ([371) allows to rewrite Eq. ^ as [lO] 

a-i_,i, = di_,u^ + duU^ - u^d^^ - u^d^S, ~ 2r]f^^daU°' /3, 

which gives 

aox — cosh (pdt(l> — sinh (pdxcf) + cosh (f>dx£, — sinh (pdt^, (68) 
o'xx = 4cosh09a;(/)/3 — 2sinh09j;^ — 2sinh</i9t0/3. (69) 

The above expressions contain time-derivatives that make the equation ^^T'^" = a second order equation in the 
time derivative. In order to obtain a well-defined evolution equation, which is first-order in time, one has to use the 
ideal hydrodynamics equation and express the time derivatives via space derivatives. Using Eqs. (|40p - (|4ip one finds 

2w sinh 20 2?i;[l -t- cosh2(?!)] 

" "2 + cosh 20' """" " 2 + cosh 20 ' 

, ,90 . , , 9C rr^ d /sinh0\ , , 

. = cosh0--smh0^=T-(^^j. (71) 

Note, that both components involve the same combination of derivatives w. This combination contains besides the 
derivative of the velocity also the derivative of the temperature. In non-relativistic hydrodynamics such terms arise 
from the heat conduction. Here their appearance is in accord with the fact that the fiuid viscosity allows to capture in 
a unified way both the effects of viscosity and the heat conduction. With viscosity the components of the stress-energy 
tensor read 

T- . 2T^sinh20- ^^ffi^A (^j^) , (72) 
^ 2-Hcosh20 dx \ T J ' ^ ' 

T- ^ [2cosh20- 1] - 2^(A)T^[l + cosh20] 9 /sinh0^ 

^ ^ ^ 2 + cosh20 dx \ T J ^ ^ 

Thus, stationary one-dimensional solutions of the GET hydrodynamics obey (cf. Eqs. (p7|) and note that partial 
derivatives becomes ordinary ones): 

These equations pass to the ideal hydrodynamics equations in the zero viscosity limit F{X) 0. Since far from x — 
the viscous contribution in these equations can be neglected due to the smallness of the derivatives, we find that the 
asymptotic values of the functions at -l-oo and —00 satisfy 

T4^sinh20+oo - T:*^sinh20_oo = Ci/2 , T^^ [2cosh20+oo - 1] = [2cosh20_oo - I] = C2 . (76) 

These conditions with v = tanh0 correspond to the Rankine-Hugoniot conditions ((62|) . In this way the shock 
characteristics Ti and Vi in Eq. (|62p describe the asymptotic values of the fields far from the region, where viscosity 
is important. In order to construct the values of the fields in the viscous region, which corresponds to a discontinuity 
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in the shock description, one has to solve the system of Eqs. ((74)) - ([75]) . Multiplying the first equation by 1 + cosh 20, 
the second by sinh 20, and subtracting the equations we find 



^4 Ci coth 0-6*2 



r = ^ . (77) 



Plugging the above relation into Eq. (|74p we find that satisfies 
2 sinh 20 [Ci coth0- C: 



_^ i^/^F{\) d I sinh< 



2 + cosh 20 dx \ [Ci coth - Ca 



1/4 



Ci . (78) 



We seek a solution to this equation that satisfies Neumann boundary conditions, i.e. vanishing of the derivatives, at 
X = —oo and x = +oo. One can rewrite the equation as 

c/0 _ r(0)[2 + cosh20][2r4(0)sinh20-Ci] 

dx~ F(A) cosh 0[2r4(0) sinh 20 + Ci/3] ' ^ ' 

with T(0) defined by Eq. ^T^j. Wc find that v = tanh0 obeys 

dv _ i'^/^Cl^^lC - v]^/*[3 - v^][v - vi][v - V2]VT^ 



dx F(A)wi/4[^(5_-t,2) _4„] 

where ( = C1/C2 and vi^2 are the roots of ({1 + 3u^) — — 0, 



(80) 



Via = ■ (81) 

Note that vi = l/Swa, cf. Eq. (|63p . The viscous counterpart of the ideal hydrodynamics shock considered in the first 
subsection corresponds to the range 1 < C < 2/-\/3 where vi^ are real and 1/3 < wi,2 < 1- 
We first consider the case v{—oo) < v{+oo) and search for the solution of Eq. ((80)) satisfying 



2-^/1^ 2+ v/4^3C^ 
t;(-c») = U2 = ^ , u(+oo) = ui = — . (82) 

Such a solution must have a positive derivative and values close to V2 in the range V2 < v < vi in some region of 
space. However, Eq. ((80)) dictates that at least close to vi^ the sign of dv/dx is determined by [w — — W2] which is 
negative at i;2 < ■y < ui. This follows by observing that the factor C(5 — w^) — 4u in the denominator of Eq. ((80)) obeys 
C(5 — w^) — 4i> = 3([v — vi][v — V2]+4:(^[l — v'^]. Thus it is not possible to construct the solutions with v{—oo) < i;(+oo). 
In contrast, the solutions with v{—oo) > i;(+oo) for which 



2 + v/4^3C^ 2- v/4^3C^ 
t;(-oo) = U2 = , v(+oo)=vi = — . (83) 

are easily constructed. Noting that (^(5 — w^) — 4w is always positive for U2 < < i"! at C > 1 we conclude from 
Eq. ((80)) that v{x) decreases monotonically from [2 + a/4 — 3C^]/3C at a; = —00 to [2 — a/4 — 3C'^]/3C at a; = +00. 
Here we use that it takes the solution of an equation of the type dv/dt = H{v) infinite "time" t to reach a zero of an 
analytic function H {(/)). The implicit form of the solution is given by 

^ r v'^/na5-v")-4v'w 

33/4CC2'/^ A/3C [C - i^?/*[3 - v'^][v' ~[2+ v/4-3C2]/3C]K - [2 - ^4 - 3(21/3^71^^' 

where the constant xq is the point at which v{x) attains its middle value 2/3C (naturally the solution is translation- 
invariant). In the following we set xq — 0. Note that v{x, Ci, C2) = v{Ci/C2, C^^x). The dependence of temperature 
on X is determined by T^{x) = C2[C — v{x)]/3v{x). One observes that the necessary condition that the resulting 
expression for is positive is satisfied by > 1 and v < 1. The ideal hydrodynamics shock is recovered in the formal 
limit of zero viscosity, F{X) — + 0. 

Thus the viscosity selects among the ideal hydrodynamics shocks only those for which the infiow velocity v{—oo) 
is larger than the outflow velocity ^(-l-oo). This has a simple meaning: the dropping of velocity corresponds to a 
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dissipative action of the viscosity (friction) , while velocity increase would qualitatively mean an accelerating action of 
the friction. The selected flows start as supersonic and end as subsonic. 

It remains to consider the condition under which the above solution has a small Knudsen number to justify neglecting 
the higher derivative terms in and T^^ . We notice that generally F{\) is of the same order as G(A) appearing in 
the expression ([2]) for lcor^ see Ell- We define the local Knudsen number by 



Kn — L 



3C[3 - v^]\[v - vi][v - V2]\VT^ 
v[C{5 - - Av] 



(85) 



where we used Eqs. ^ and ([50]) . The viscous shock solution above is valid if Kn ^ 1 uniformly. At fixed ^, 
by choosing a sufficiently small C2 one may achieve an arbitrarily large scale of variations of v in space, due to 
v{x, Ci, C'2) = v{C'i/C2, C^^^x). This, however, does not lead to a solution with a small Knudsen number, since the 
correlation length also increases under such scaling of d . This feature is taken into account by Eq. ([85]) , that depends 
only on the ratio Ci /C2 . 

To check that Kn ^ 1, it is enough to substitute in the above equation the middle value 2/3^ of v which produces 
a term of order 4 — 3^^ {( belongs to the range [1, 2/-\/3]. It follows that Kn <C 1 is satisfied only for C close to 2/^/3 
(a more precise condition is derived below). For such (, the inflow velocity w(— 00) is close to the speed of sound 
l/VS, and the outflow obeys ?;(+oo) « v{—oo), ci. remarks to Eq. (|63p . The Knudsen number is small just because 
the total variation of v is small. In other words, for this solution the requirement of locally small variations of fields, 
implies also that the global variation is small. 



D. The viscous shock solution at a small Knudsen number 



In the small Knudsen number domain — 2/\/3 — e with e ^ 1, the velocity changes in the interval [l/-\/3 — 
(5, l/\/3 + S] where S = e^/^3^^/''. Using that the total variation of the solution is small, one can obtain more explicit 
expressions in this case. To leading order in 5 Eq. (|80p becomes 



dx ~ F{\) 



(86) 



where v = v — l/\/3 is the velocity deviation from the speed of sound. The solution of the above equation gives the 
well-known kink solution, which can also be obtained by expanding the integral 



V = — = — 5 tanh 
V3 



L = 



F{\) 



31/4V2C2 

where we used that T « [6*2/3]^/'* holds for the solution. One has 

C2 , 2C2(5 



1/4, 



cor 



V3 



tanh 



(87) 



(88) 



We see that the Knudsen number Kn — Icor/L ~ (5 is small at 5 ^ 1. Note, that the Knudsen number definition ((85|) 
would give incorrectly Kn ^ - the reason is that v{x) has a special form where v'{x) is determined by the small 
correction v to the main constant part 1/V3- 

One can exclude the constant component in v(x) by passing to the frame moving at the speed of sound, where the 
solution takes the form of a kink propagating at the sound speed. 



tanh 

2 




ji4 ^2 



C2 

3 



2C2S 
V3 



tanh 




(89) 



in the leading order in S, where we restored c in v. One sees that the solution is a near equilibrium state. In the limit 
c — *■ 00 the motion gets degenerated to v = const in accord with the fact that the only incompressible one-dimensional 
flow is a constant one. 
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E. The shock solution at a Knudsen number of order unity 

We saw above that if the shock solution is such that the inflow velocity w(— oo) differs significantly from the speed 
of sound then there are regions in the flow where Kn 1 (this follows from self-consistency: if Kn <C 1 then viscous 
hydrodynamics should work well) . The solution of the viscous hydrodynamics then does not describe the solutions 
correctly in those regions though it may give a good idea of the solution there. To find an exact solution one should 
account for higher order derivatives in the constitutive relation for T^"^ (for instance, the second derivative terms were 
discussed in [lB|)- The above however does not undermine the utility of the ideal hydrodynamics shock which is 
still a good description of the solution far from the transition region near x = 0. 

We show this point explicitly. We consider a non-equilibrium steady state of a CFT fluid which is sustained by the 
fluxes of energy and a;— component of momentum from infinity. In this state the mean value of T^" depends only on 
X and Eqs. (|67|) are satisfied irrespective of the validity of the hydrodynamic constitutive relations. Now, far from the 
transition region the scale of variations of the mean value of the stress-energy tensor becomes large, the state is locally 
very close to equilibrium and the ideal hydrodynamics approximation T^^, = [rj^v + ^u^u^,] works with increasing 
degree of accuracy. In other words, far from the transition region the state looks locally closer and closer to the 
global equilibrium (the information that the local equilibrium state is not the global one becomes further separated in 
space). As a result, the ideal hydrodynamics shock solution allows to establish correctly at all orders the connection 
between the asymptotic values of w(— oo), r(— oo) and v(-l-oo), T(-I-cxd), and thus also of T^^'^{—oo) and T^^'^{+oo). As 
far as the intermediate region of the transition is concerned, it is described by the solution in the previous subsection 
for the inflow velocity close to the speed of sound l/\/3. When the inflow velocity l/VS < oo) < 1 is not close to 
l/\/3 the description of the transition region demands the account of the higher order derivatives. We do not expect, 
however, that the higher derivative terms will lead to a further selection of the ideal hydrodynamic shocks beyond 
the one imposed by the viscosity. 

F. How is it possible to have Kn <^ 1, while the viscous term still being important? 

It is natural to ask why is it possible to construct a solution where the viscous term, whose effect is generally 
measured by the Knudsen number, is important while higher order derivative terms, whose effect is measured by the 
same parameter, are not. The reason being the presence of an additional dimensionless parameter in the equations, 
v/c (here we restore c in the relations). As long as we consider relativistic motions with u ~ c, the viscous to ideal 
terms ratio is governed by the same parameter Kn as the ratio of higher order derivative terms to the viscous one. 
However, as we pass to motions with v c the viscous to ideal terms ratio involves an additional small parameter 
V /c. In the case of properly defined non-relativistic motions obeying the incompressible Navier-Stokes equations the 
dimensionless parameter governing the viscous to ideal terms ratio is the inverse of the Reynolds number Re 

Re - vL/iy , (90) 

where v is the characteristic velocity, L is the characteristic scale and i' is the kinematic viscosity, sec Eq. (IT^ . Thus, 
one can have a situation, where in some regions of space the Reynolds number is not large and viscosity is important, 
while the Knudsen number is still small everywhere. As long as the viscosity prevents the generation of smaller scales 
by the dynamics (as it does at least in known cases), one obtains a global solution of the hydrodynamics with viscous 
term important but higher order derivatives not. In other words, while viscosity is typically a singular perturbation 
producing order one effects in the solution [l7j . the higher order terms are often regular perturbations. 

The relevance of the above discussion to the viscous shock solution considered before, is that at small Knudsen 
number, v{x) <^ c in (|89p . Thus, in the frame moving at the speed of sound one deals with the situation of small v/c 
and the above picture holds. Let us stress again, however, that this solution has no well-defined non-relativistic limit 
as the shock propagates at a speed of order c. 

Note also, that here and below we use the notion of a local Reynolds (Knudsen) number, which is the Reynolds 
(Knudsen) number defined via the characteristic values of the scale and velocity in the considered region of space . 

V. GRAVITY AND FLUID DYNAMICS 

It has been shown in 3], that the four-dimensional CFT hydrodynamics equations are the same as the equations 
describing the evolution of long-scale perturbations of the five-dimensional black branc. In fact, the derivation of this 
result parallels exactly the conventional derivation of hydrodynamics equations from the Boltzmann equation. Let us 
show this in the simplest case of a non-relativistic monoatomic gas. 
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The Boltzmann equation for a non-relativistic monoatomic gas reads [18| 

dtf + VW^f^Stf, (91) 

where f{x^p,t) is the density function in a single particle phase space, V — p/m where m is the particles mass, 
and Stf is the collision integral, which is quadratic in /. Equation ([OTjl has a steady-state solution of a Maxwell 
distribution form 



[zTTmkB-L \ ' 



m {V — v)^ 
2kBT 



(92) 



where N is the mean concentration of particles, T is the temperature and v is an arbitrary velocity present due to 
the Galilean invariance of Eq. (pij) . 

The hydrodynamics equations arise when solving the Boltzmann equation by the method of variation of constants. 
Consider a solution in the form 



N{x t) 

f{x,p,t) = fo{x,p,t) + 6f{x,p,t), fQ{x,p,t) = ^2TrmkBT{x t)]^/^ 



2kBT{x,t) 



(93) 



where 5f <C /o- This describes the local thermal equilibrium in the gas. The construction of a solution to the 
Boltzmann equation where /o is the zeroth order approximation, produces a series in the Knudsen number Kn^ i.e. 
the ratio of the mean free path to the scale of spatial variations of the fields. Thus, 5 J ^ Knfo, and the hydrodynamics 
equations governing the evolution of N{x,t), T{x,t) and v{x,t) result as the conditions of the series constructibility. 

Analogous derivation of hydrodynamics from the (relativistic) Boltzmann equation can be performed for a CFT. 
This becomes possible in the limit of the weak coupling where quasiparticle excitations are well-defined (see e.g. 0]). 
In contrast, at a strong or order unity coupling the Boltzmann equation does not apply. What this means normally is, 
that for non weakly coupled fluids the hydrodynamics equations are valid but cannot be derived from a microscopic 
description. The CFT fluids, however, bring a surprise. 

Consider the derivation of the CFT hydrodynamics at strong coupling from the equations of gravity. The analogue 
of Eq. (PT|) is given by the five-dimensional Einstein equations with negative cosmological constant 

+ 4<?„„, = 0, R=-20. (94) 

These equations have a particular "thermal equilibrium" solution - the boosted black brane 

ds^ = -2Uf,dx^'dr + n'^T'^r-^Uf.u^dxf'dx'' + r'^'q^.^dx^'dx" , (95) 

where the temperature T is constant. The presence of an arbitrary constant vector in the above solution is due to 
Lorentz invariance, much like the presence of v in Eq. ([5^ due to Galilean invariance. As in Eq. ([M]) . one looks for 
a solution of the Einstein equation by the method of variation of constants using the ansatz 

9mn = (.go)mn + Sg^n, (so)™™^?/™^?/" = ~2u ^{x'')dx^' dr + n^T\x")r~^u^{x")u,{x")dx''dx'' + r^r^^.dx^dx'^ ,{96) 

where y = {x^,r). As in the Boltzmann equation, the condition of constructibility of the series solution produces 
equations for u'^{x") and T{x°'). These equations are the hydrodynamics equations for a certain class of CFTs in the 
limit of infinite coupling Q . As far as the hydrodynamic equations are concerned the infinite coupling means only a 
particular value of constants in the relations. For example, for dissipative hydrodynamics considered in the previous 
sections infinite coupling means taking instead of F{X) its appropriate A ^ oo limit, F = I/tt. Note, that at the level 
of the ideal hydrodynamic equations there is no difference at all between different CFTs - the equations are the same. 
The resulting series for gmn is the series in the Knudsen number of the boundary CFT hydrodynamics. 

Thus equations of gravity play for the hydrodynamic equations of a strongly coupled CFTs much the same role 
as the Boltzmann equation plays at a weak coupling. It is natural to inquire whether the equations of the CFT 
hydrodynamics allow for a microscopic derivation also at an arbitrary intermediate coupling (say from the equations 
of gravity with higher order curvature terms in the action). This important question is beyond the scope of this work 
and is a subject for the future research. 

We conclude from the above that the solutions to the CFT hydrodynamics at an infinite coupling describe not 
only the large-scale motions of the theory, but also the solutions to gravity corresponding to Eq. ([96]) . We may 
now use the general results for the CFT hydrodynamics, obtained in the previous section and valid in particular for 
the hydrodynamics of the infinitely coupled CFT resulting from Eqs. (|94p. in order to gain some insight into gravity 
dynamics. If we consider initial conditions for gravity, that correspond to the initial conditions for hydrodynamics with 
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V c and Kn <C 1 everywhere, then for some Umited amount of time we may use ideal hydrodynamics to determine 
the metric evolution. The ideal hydrodynamic flow is compressible, and it leads to the formation of discontinuities 
(and maybe other, yet unknown types of singularities) in a finite time. In particular, for "one-dimensional" initial 
conditions, one may use the exact solutions obtained in the previous sections. Then, ideal hydrodynamics predicts 
that in a finite time, regions with large metric curvature will appear. They correspond to large values of gradients 
of the hydrodynamics fields. Eventually, in these regions higher order derivative terms will become important for the 
evolution of and T and (expectedly) stop the growth of the curvature. One expects Kn 1 and Sg ^ go in these 
regions. 



A. The importance of the solutions to non-relativistic hydrodynamics for the construction of global 

solutions to gravity 

As already explained in the last subsection of section IIVI the consideration of non-relativistic solutions with f <C c 
allows one to construct global solutions of hydrodynamics with Kn <ti 1 everywhere. Such solutions are especially 
important from the viewpoint of the construction of globally valid approximations to the solutions of gravity. Indeed, 
for them one has the explicit form g w go of the approximate solution to the Einstein equations (|94p . This is quite 
similar to the situation for non-relativistic classical gas where the local thermal equilibrium / ~fo holds even though 
viscosity is important for the evolution of the local parameters of equilibrium T, N and v, cf. [la|. 

The explicit form of the approximate solution to the Einstein equations resulting from the solution to the non- 
relativistic incompressible Navier-Stokes equations was summarized in Eqs. (|17p -(|18 p . 



B. Global solution for gravity from a small Knudsen number viscous shock 



As an example of a global solution to gravity obtained from hydrodynamics we consider the five-dimensional 
asymptotically AdS gravity solution that corresponds to the viscous shock. Using the results of the previous section 
we find that the metric g^s defined by 



igvs)mndy'^dy"- = V6 



\/3S tanh(a::/L) 



dtdr - V2 



3V3(5 tanh(a;/_L) 



dxdr + r^rj^iydx^dx'^ 



3r^ 



1 + tanh 



1 — 3\/3(5 tanh 



L 



(!) 

dx^ 



L 



\/i5 tanh ( — 



- V3 1 - 2\/35tanh 



(!) 



dxdt 



(97) 



is a solution of the Einstein equations to order 5^ <C 1 (the constant C2 > is arbitrary), and indeed, an explicit 
calculation verifies that. 



Above we used that at small S one has uq — — [I — w^] 



21-1/2 



V35tanh(a;/L)/2] and 



v[l- 



^1/2 



l/V2[l - 3V35tanh(a;/L)/2], where v{x) is given by Eq. (|87p . As we saw in the previous section the requirement 
Kn <^ 1 imposes that the total variation of v is small and the above solution describes a near equilibrium state. 

The above stationary solution corres pon ds to the solution of the Einstein equations with the boundary condition 
imposed on the boundary stress tensor [10| 



Ti: - -2 hm {K^ - 6^^) 



(98) 



where K^i, is the extrinsic curvature tensor to the surface at fixed r and T^'^ is the solution of the hydrodynamics 
equations. If one used the ideal hydrodynamics approximation, where the shock is discontinuous one would obtain 
a domain wall solution for gravity with no 5— function sources at the wall. At small Kn, the solution l|97[) allows 
to resolve the inner structure of the wall. At Kn 1, the resolution of the wall structure requires the inclusion 
of the higher order derivatives in the analysis. However, we stress again that the ideal shock solution provides a 
correct description of the domain wall solution far from wall at any Kn. Finally, we note that it is straightforward to 
generalize this domain wall solution of gravity by applying special conformal transformations to the solution. 
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